Transport spectroscopy of NS nanowire junctions with Majorana fermions 
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We investigate transport though normal-superconductor nanowire junctions in the presence of 
spin-orbit coupling and magnetic field. As the Zeeman field crosses the critical bulk value Be of 
the topological transition, a Majorana bound state (MBS) is formed, giving rise to a sharp zero- 
bias anomaly (ZBA) in the tunneling differential conductance. We identify novel features beyond 
this picture in wires with inhomogeneous depletion, like the appearance of two MBSs inside a long 
depleted region for B < Be. The resulting ZBA is in most cases weakly split and may coexist with 
Andreev bound states near zero energy. The ZBA may appear without evidence of a topological gap 
closing. This latter aspect is more evident in the multiband case and stems from a smooth pinch-off 
barrier. Most of these features are in qualitative agreement with recent experiments [Mourik et al, 
Science, 336, 1003 (2012)]. We also discuss the rich phenomenology of the problem in other regimes 
which remain experimentally unexplored. 



Following early ideas based on exotic p-wave super- 
conductors m [2:, it has been recently predicted that 
Majorana quasiparticles should appear in topological in- 
sulators |3j and semiconductors with strong spin-orbit 
(SO) coupling [IHS]. In proximity to s-wave superconduc- 
tors, these systems behave as topological superconduc- 
tors (TS) when the excitation gap is closed and reopened 
again: as the gap crosses zero, Majorana bound states 
(MBSs) appear wherever the system interfaces with a 
non-topological insulator (see Refs. [9^ 10 for reviews). 

The TS transition occurs when an external Zeeman 
field B exceeds a critical value Br 




fined in terms of the Fermi energy /i and the induced 
s-wave pairing A [7(181. This prediction has spurred 
a great deal of experimental activity towards detecting 
MBSs in hybrid superconductor-semiconductor systems. 
Indeed, signatures of Majorana detection have been re- 
cently reported in Ref. [11]. These experiments (and 
subsequently Refs. [l?, T3^) clearly show the emergence 
of a zero-bias anomaly (ZBA) in differential conductance 
dl/dV measurements as B increases. It has been pre- 
dicted that such ZBA reflects tunneling into the MBS 
[MHTBI. Crucially, the emerging ZBA, which signals the 
TS transition, should be accompanied by a closing and 
reopening of the excitation gap [TT, TB , something which 
is however not observed. Other experimental findings, 
such as ZBA splitting and coexistence of Andreev bound 
states (ABSs) and MBSs [11 , need also further analysis. 

Motivated by this, we present here a detailed study 
of transport through normal-superconductor (NS) junc- 
tions containing topological wires. As in the experiment, 
the wires are tunnel-coupled to the normal reservoir by a 
pinch-off gate Vp (to allow for transport spectroscopy us- 
ing dl/dV)^ and are depleted by an additional gate Vd in 
order to bring Be down to accessible fields (since the in- 
duced potential Ud lowers the Fermi energy ji ^ (i — Ud). 
The depletion profile, however, is necessarily inhomoge- 
neous^ since it cannot extend deep into the supercon- 
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FIG. 1. (Color online) Schematics of the nanowire junction 
in the N'^S'^S (a) and NS"^ S (b) setups, and spatial varia- 
tion of superconducting gap and potential profiles (c). Gate 
Vd depletes the wire, while Vp creates a tunnel contact (I) to 
the left (normal) reservoir. One (red) or two (red and yellow 
spheres) Majorana bound states may appear at the edges of 
the depleted region depending on the Zeeman field and gate 
voltage Vd. (d) Transport regimes for a transparent NS junc- 
tion (Vd,p = 0, /i = 4 A) in the Zeeman field - bias plane. 



ducting side due to efficient screening, see Fig. [T] and 
Ref. |TT| We find that for this class of devices, a number 
of distinct transport regimes arise as the various sections 
of the wire transition to different electronic phases. We 
characterize these regimes and the rich phenomenology 
that results beyond the simplest picture [MHTB] . In par- 
ticular, we address the question of whether the ZBAs 
are related to Majorana physics, and confirm that this 
is indeed the case for long depletion regions. We also 
analyze the development of ABSs close to zero energy 
when the pinch-off gate lies at a finite distance from the 
NS junction. Our main results are summarized in Fig. 



^e) where we demonstrate that the dl/dV of reahstic 
junctions with inhomogeneous depletion and multisub- 
band fihing may not show a distinct closing of the gap 
and yet exhibit ZBAs of Major ana origin. In most cases, 
these ZBAs show a residual splitting and may coexist 
with ABSs, features also observed in Ref. 11, 

Model. — We first consider a one-dimensional NS junc- 
tion [see Fig. Ilia)], with a BCS-type Hamiltonian H = 



Ho^H, 



pairing? 



where 



^0 



Jd.^Hx) 



2m 



iacjydx + B(Jz + U{x) — ji 



iIj{x) 



and impairing = J dx ilj\x)iA{x)ayilj\x) + h.c. 

Here a is the SO coupling and B is the Zeeman split- 
ting (given hy B = gjaB^/^^ where B is an in-plane 
magnetic field, jub is the Bohr magneton and g is the 
nanowire g- factor). We assume a position dependent 
pairing A{x) induced by the superconducting electrode 
such that A{x ^ oo) = A and A{x -^ -oo) = 0. The 
term U{x) = Ud{x) + Up{x) is composed of two parts: 
Upix) comes from the pinch-off gate Vp in the normal re- 
gion at a distance L^ from the NS interface, and Ud{x) 
models the potential induced by the depletion gate Vd 
[19] . Gate Vd may extend all the way into the nor- 
mal side of the NS interface [case N^S^S, with a de- 
pleted length Ld = L^d -\- Lgd^ Fig. flla)], or be lim- 
ited to the end of the superconducting side [case NS^S, 
Ld = Lgd^ Fig. [ijb)]. We will consider the former case 
first, where we cover different parametric regimes, and 
then turn to the second one, which is closer to the ex- 
perimental setup [11 . Realistic experimental parameters 
are: A = 250/ieV is the induced gap that, for an InSb ef- 
fective mass m = O.OlSmg, corresponds to a length scale 
I/A = h/VmA = 142nm. Strong SO coupling, represen- 
tative of InSb wires [20], is a = 20 meV nm, with SO 
length Lso = h^/{ma) = 200nm = 1.4La pT]. 

Scales. — A localized MBS is formed at the boundary 
of a trivially gapped and a TS portion of the wire. At a 
point X the wire will be in the TS phase if A{x) > and 



B> V[/i-/7(x)]2 + A(x)2 



(1) 



The asymptotic value of the critical field is the proper 
(bulk) critical field Be. Apart from Bc^ several other 
Zeeman scales dictate the junction's transport properties. 
The first one is the TS critical field in the depleted part of 
the superconducting wire, B^ = ^/{Jl^^TJdf^^r~^ .> that 
is smaller than 5c, as is the purpose of the depletion gate. 
It should be noted, however, that the depleted S^ re- 
gion has a finite length, which crucially affects Majorana 
modes for B^ < B < Bc^ diS discussed later, while the S 
portion is assumed infinite. Secondly, there is the field 
above which the normal side of the wire becomes a helical 
liquid (momentum and spin become correlated). In the 
NS^S case (normal side not depleted), this is B^ = /i, 
which is typically slightly smaller than Bc^ but bigger 
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FIG. 2. (Golor online) Density plots of the dl/dV in the 
N'^S'^S junction (/x = 4A, Ud = 3.25A, Up = 25A, ^ = 0) 
for Lso = 1.4La as a function of bias voltage V and Zeeman 
field B with a tunnel pinch-off barrier and a depletion region 
of length L^d + L^d, Fig. [ija). Different columns feature 
increasing values of L^d from left to right, whereas different 
rows feature increasing length Lj^^d from top to bottom. 



than both B^ and the corresponding helical field in the 
N^S'^S case, namely Bf^ = \fi-Ud\ < B^. Finally, there 
is the superconducting gap itself, B/^ = A, whose signif- 
icance will become clear later. All these scales {Ba plus 
Bf < Bh < Be in the NS"^ S case, or Bf^ < B^ < Be in 
the N^S^S)^ control different aspects of the junction's 
differential conductance in the B-V plane. 

Differential conductance. — The dl/dV of a NS junc- 
tion may be related to the intrinsic conductance at zero 
temperature by the expression [22] : 



dliV) 
dV 



e 
~h 



Ar-Tr(rLree) + Tr(4re/,) 



Here, M is the number of propagating channels in the 
normal side at energy e = F, and Vee and Veh are their 
normal and Andreev reflection matrices. These matrices 
can be computed in a number of ways. The most flex- 
ible is the recursive Nambu Green's function approach, 
employed here (for full details, see [23]). 

Before considering the effect of t/(x), we show the 
transport phase diagram [see Fig. [l]^c)] in the simple 
NS transparent limit, i.e. in a regime where the concepts 
of MBSs and ZBAs do no longer hold. We observe dif- 
ferent transport regions in the B-V plane characterized 
by an integer dl/dV ^ ne'^/h^ with n = 0, 1, 2, 3, 4. Such 
is the case of Cooper pair transport (region I, n = 4) or 
single quasiparticle transport (region III, n = 2). The 
latter is a TS regime, whose topology becomes evident 
in the dl/dV despite the fact that the associated Majo- 
rana fermion is completely smeared out due to the gapless 
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FIG. 3. (Color online) dl/dV in a short (a) and long (b) 
N'^S'^S junction (/x = 4A, Ud = 3.25A, Up = 22A, ^ = 0). 
Coloured guidelines represent the gap at large (green) and 
small (blue) momentum (see bandstructure in the left inset) 
in the depleted (dashed) and bulk (solid) superconducting re- 
gions. Vertical solid (dashed) line marks Be (Be), (a) For 
short depletion region, a mini gap transport regime is devel- 
oped for B > Ba (green dot), unrelated to Majorana modes, 
(b) Long depletion region, however, allows for the existence 
of ABSs around zero energy below B^ (red dot) and (split) 
Majorana fermions for B < Be, see blowup in the right inset. 



spectrum for x < [24H27] . Between these two regions, 
the helical regime is characterized by a fully suppressed 
zero bias conductance (region II, n = 0). These results 
extend the concept of half-integer conductance quantiza- 
tion [24] beyond linear response. 

We now consider the N^S^S junction with the full 
U{x). Its dl/dV response (with Lso = 1.4La) is plotted 
in Fig. |2j Different panels cover different ratios L^d/L/\ 
and Lsd/L/\- The tunnel barrier Up is tuned in each 
case to yield spectroscopic resolution in the transport 
response. A wide range of behaviours become appar- 
ent, that reflect the local density of states (DOS) at the 
pinch-off gate. The most paradigmatic one is probably 
the one in the top-left panel. It reflects the closing of the 
effective superconducting gap (marked by the gap-edge 
conductance peaks) as B increases. The gap-edge DOS 
peak transforms into a Majorana mode at zero energy for 
B > Be [pictured in Fig. flla) as a red sphere]. The gap 
reopens in the TS phase (see solid blue line in Fig. [3|, but 
the local DOS at the contact is no longer peaked because 
the spectral weight is transferred to the Majorana mode. 
Hence, a prototypical three-pronged structure arises in 
the B-V plane. However, the relevant phenomenology 
is by no means exhausted by this. Different scenarios 
arise at large Lj^d (Fig. ^ bottom row), with the devel- 
opment of ABSs, or large L^d (right column), with the 
development of ZBAs below Be. Further phenomenology 
is obtained by varying Lso [23] , 

Short N^S^S junction. — We now analyze in more de- 
tail the short junction case, Ljyd^Lsd <C La (top-left 
panel in Fig. ^. The scale B^ has little meaning in 
this case, since any pair of Majorana modes forming at 
the ends of a short TS wire will strongly hybridize into 



two conventional states with energies ~ ±A. The effect 
of increasing Lso is to flatten the gap-edge dl/dV peak 
(bright yellow) to lower energies, as shown in Fig. IsFa) 
(Lso = 5La). At these large Lso^ spin becomes a good 
quantum number and the Zeeman field splits the particle 
and hole bands of the S region into spin up/down sub- 
bands (white arrows in Fig. |3|. Particle and hole gap 
edges anticross at zero energy for B > B^ (green dot), 
resulting in a minigap (which vanishes for Lso ^ cxo) 
with almost flat edges near zero. This structure appears 
as a (split) ZBA below Be that is unrelated to Majorana 
mode formation [28] . 

Long N^S^S junction. — Increasing the length Ljyd of 
the depleted normal wire gives rise to ABSs in the N^ 
region that are probed by the tunnel barrier Up (see bot- 
tom row of Fig. [2]). Increasing Lso we find these ABSs 
resonances approaching a degenerate zero energy cross- 
ing, see Fig. [SJ^b). These low energy ABSs, however, 
cease after the helical transition at Bf^ (red dot). While 
the ABSs move up and down in energy for B < Bf^ de- 
pending on their spin character, as soon as the N^ region 
becomes helical, all ABSs disperse away from zero energy 
with increasing 5, since incoming states do not have spin 
partners anymore. Beyond the second threshold Bf (ver- 
tical dashed white line), the local band gap closes in the 
gd ^Qg-^QYi (dashed blue line) and a ZBA builds up, which 
is caused by Majorana fermion pairs forming in the S^ 
region of length L^d, also assumed large. The finite L^d, 
however, produces a residual splitting of the two Majo- 
ranas, visible in the ZBA that oscillates with B [29H32] 
as long diS B < Be. Above the bulk Be, the rightmost 
Majorana escapes to x ^ 00 (see Fig. [l]), where it no 
longer overlaps with the one at the A^^S*^ interface, and 
the splitting vanishes [see inset of Fig. [SJ^b)]. 

NS^S junction. — We now turn to the type of setup 
explored in Ref. [11] the NS^S case. The crucial differ- 
ence with the N^S^S setup is that the ABSs zero energy 
anticrossings may coexist with Majorana ZBAs in the 
gd region, since the N region becomes helical after the 
gd YQgiQYi becomes topological. When that happens, the 
zero energy ABSs repel the MBS wavefunction (usually 
delocalized into the N region) back into the S^ region, 
hence decoupling it from the lead. As a result, the ZBA, 
as measured by the dl/dV, is suppressed (see arrows in 
Fig. [4J where Lso = 1.4La). In a single mode wire (Fig. 
pj top row), this may happen only for B < Bh < Be, 
since it requires an interface between a non-helical N and 
a topological S^ [33]. Interestingly, the above behavior 
agrees with the experimental observation of an intermit- 
tent ZBA that disappears and reappears as an ABS anti- 
crosses at zero energy. Constant-5 sweeps of Ud are also 
found to closely correlate with the experiment [23] . 

A further feature apparent in the experiments is the 
absence of gap-edge singularities closing just before the 
formation of the ZBA. Up till now, all potential pro- 
files have been assumed spatially abrupt [decay length 



(5 = in all profiles of Fig. flic)]. Imperfect screening, 
however, will lead to gate-induced potentials that decay 
slowly along the wire, specially at low electron densities. 
When a smooth Up{x) profile is taken into account, the 
gap-edge peaks are quickly washed out, and the dl/dV in 
the tunneling regime is no longer a perfect measure of the 
local density of states [34] . Unlike for sharp pinch-off bar- 
riers, transmission through smooth barriers mostly pre- 
serves longitudinal momentum. Since moreover barrier 
transmission is lower for smaller momenta, the gap-edge 
resonance closing around zero momentum at B^ is poorly 
probed by a smooth barrier. Visibility is restored as B 
increases beyond B^^ since then the Zeeman splitting in- 
creases the momentum components of all resonances, in- 
cluding the ZBA. Fig. [4lb) shows this effect. Note on 
the other hand that the gap-edge signal at higher ener- 
gies ~ A is roughly constant in B. This corresponds to 
the large-momentum band edge, not the true gap edge at 
small momentum. The former band edge shows no sign 
of the different transitions (the TS transition at Be in 
particular), and never closes. Therefore, all zero energy 
structure appears disconnected from any closing of the 
gap, as measured by the dl/dV. 

Transport features connected to large momenta should 
be enhanced in a multisubband systems. The pinch-off 
condition for tunneling spectroscopy requires a higher Up 
barrier in this case, to shut off the additional open modes, 
which will necessarily contribute with a stronger signal 
relative to the shallower, lower momentum mode. This 
is likely the experimental situation in Ref. [11] where 
the observed gap edge signal exceeds the amplitude of 
the ZBA. In order to confirm this idea, we have per- 
formed multisubband transport simulations (see imple- 
mentation details in [23i). In Fig. [4jc), (d), we show 
results corresponding to a fixed /i defined from the bot- 
tom of the topmost subband, which becomes topological 
at low magnetic fields, while the lower ones remain triv- 
ial. The single mode analysis carries over unchanged, 
albeit with a stronger gap edge and the associated high- 
momentum structure of the additional modes. Fig. ^e) 
shows constant B traces of Fig. ^d) at finite tempera- 
ture {ksT = 70mK). Thermal smearing washes out the 
ZBA splitting, whilst also reducing its height. 

Conclusions. — We have identified various transport 
regimes in depleted NS nanowire junctions with SO cou- 
pling and Zeeman field. Depending on the Zeeman field, 
the wire as a whole may be in any given combination 
of helical/non-helical and trivial/topological phases for 
its normal and superconducting portions, respectively. 
These include a helical depletion {Bf^ < B < B^) and a 
topological depletion {B^ < B < B^^Bc) phases that 
arise at fields below those of the proper helical bulk 
{Bh < B < Be) and topological bulk {B > Be) regimes. 
The different phases have distinct sub-gap signatures in 
transport, particularly if the depleted S^ region is long 
enough. In this case, ZBAs appear that are caused by the 
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FIG. 4. (Color online) dl/dV response at zero temperature 
of an NS^ S junction similar to that of Ref. 11, for one open 
(spinful) mode (a,b) and two modes (c,d). ABSs coexist and 
suppress (white arrows) the ZBA from Majorana pairs within 
the depletion region. For increasing smoothness of the pinch- 
off profile Up{x) (b and d), the local gap edge peak is washed 
out, leaving an isolated ZBA ai B < Be without a visible 
closing of a transport gap. Vertical guidelines and color scale 
like in Fig. p] (e) Finite temperature dl/dV version of (d). 
Different curves are shifted for clarity from B — (bottom) to 
B = O.GBc in 40 steps. Parameters: /j. = Ud = 8A, Up ^ 20- 
30A, Ln = 450nm, L^ = l/^m, 6 = lOOnm, B^ = O.UBc 



formation of either a single (in the topological bulk phase) 
or a pair {topologial depletion phase) of Majorana modes 
in the junction. The latter is characterized by a residual 
ZBA splitting. In the case of short junctions, Majoranas 
cannot develop below Bc^ although one may still distin- 
guish between the conventional and minigap transport 
regimes (unrelated to Majorana physics) if Lso is large. 
Apart from this general analysis, we have discussed 
a configuration similar to the one in the experiment of 
Ref. [TTl We argue that at least a number of non-trivial 
features observed in this experiment are consistent with 
most of our results corresponding to transport through a 
multi-mode non-helical normal/topological depleted su- 
perconductor/trivial bulk superconductor junction, host- 
ing Majorana fermion pairs within the central region 
[Fig. p[d), (e)]. While a qualitative correspondence can 
be traced between our observations and the experiment, 
a more quantitative agreement is beyond the goal of our 
study. This is in part due to a considerable number of 
unknown parameters (precise Lso^ screened gate poten- 
tial profiles, localization and pair-breaking effects in the 
bulk superconductor [35], etc.). Moreover, other features 
that we find, such as the ZBA splitting oscillations in B, 
have not been observed. Thus, a different physical ori- 
gin of the measured ZBAs (reflectionless tunneling [36] , 
Kondo, etc) cannot be completely excluded. The lat- 
ter, however, should be accompanied (at least in its most 



conventional form) with even-odd effects as a function 
of gate voltages ^7\ , which are not apparent in Ref. [11] 
Disorder (not considered here) might also be of relevance. 
However, its precise role on Majorana physics is currently 
under active investigation [38-40J. 
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I. NUMERICAL METHOD: RECURSIVE 
GREEN'S FUNCTION ALGORITHM 

The Hamiltonian of the system with a single mode is 



H = Hq 



^pairing, where 



^0 



dx '0' (x) 



2m 



- iacFydx + BcFz + U{x) — ji 



tlj{x) 



and impairing = J dx ilj\x)iA{x) ayi/j^ {x) -\-h..c are defined 
on an unbounded one dimensional spatial continuum, 
times an SU(2) subspace for spin. A strategy to solve 
elastic transport is to discretize the continuum into an 
infinite tight-binding chain with lattice spacing a, where 
derivative dxi^ix) is approximated by finite differences. 
The computed transport properties should be indepen- 
dent on the discretization parameter, and should there- 
fore be constrained to low energies. 

The way to deal with the pairing potential is to in- 
troduce Nambu space notation. Particle and holes are 
treated much like spin, as an extra quantum number. 
This is valid in the absence of additional interactions, 
and introduces a doubling of the Hilbert space. Wave- 
functions and propagators are consequently doubled in 
size. Other than that, the recursive Green function al- 
gorithm to computing the scattering matrix^ is the same 
than in systems without superconductivity. We summa- 
rize it in what follows. 

The scattering matrix is computed by solving the 
Lippmann-Schwinger equation^, that in the form most 
useful to us relates the full scattering states in the wire 
l^a) at energy uj to the scattering states of the semi- 
infinite leads 1^^), the Green function G{uj) of the full 
wire, and the operator V coupling the leads to the sys- 
tem. 

|*«) = |*°) + (GFWie) 

The separation of the scattering problem into an easy 
lead subproblem, coupled to a complicated system is 
common. It requires identifying a repeating block in the 
Hamiltonian at x ^ ±oo (a unit cell), each of which has 
matrix elements h under H^ and are coupled to adjacent 
blocks by matrix elements v. Everything else falls into 
the system portion, and transport through it is computed 
recursively with the boundary conditions imposed by the 
leads. 

To solve the leads, one assumes they are semiinfinite 
and decoupled from the system, and uses any method 
available to solve their Green function g for a given h 



and V, evaluated on the last unit cell, the one coupled to 
the system. This can be done by self-consistent iteration 
(imposing that adding a new unit cell to a semi-infinite 
lead does not change the surface Green function g) , or by 
turning the problem into a generalized eigenvalue prob- 
lem. In any case, g should be obtained under the bound- 
ary condition that it be non-divergent as one advances 
far into the leads (it should be the retarded g). 

The next step is to slice the discretized system portion 
into subparts analogous to the unit cells of the leads. 
Unlike for the leads, these pieces will not be identical, 
and there will be a finite number of them. They should, 
however, be coupled only to nearest neighboring slices. 
It is always possible to slice the system in this way. To 
return to the scattering matrix, one needs the full G of 
the leads+system problem. This is obtained by including 
in the system Hamiltonian the self-energy E = VgV^ 
induced on the surface slices of the system by the two 
(or more) leads, which reduces the problem to the system 
alone (leads are incorporated exactly through S). This 
suggests a good way to define 'slice zero' of the system, 
i.e. the surface slice coupled to the leads. 

Once the leads are integrated out, the full propagator 
Goo within the surface slice of the system is computed 
recursively. Step zero is a simple matrix inversion 



G. 



(0) 
00 



(u; - iJoo - S)-\ 



where H^ra are the matrix elements of H between slice n 

(k) 

and m. Likewise Gnm is the propagator between slice n 
and m including the k first slices. The next step is to add 



the next slice through the Dyson equation, Gqq 



(1) 



G, 



G^'n^i/0l(c^ 



(0) 
00 " 



.00 -01V- - Hii)-^HioG^^J. The term H{^ - H)-^H 
is the self energy brought in by the additional slice. This 
can be repeated until all slices are exhausted. The general 
recursive relation is 

^00 -^00 ^^o^k-l-'^k-l,k[(^- -tikk) J^k,k-l^k-l,0 

The recursion requires to keep track of Goo, Go/c and 
G/eo in each step. The computational cost scales linearly 
in the number of recursion steps (the number of slices) 
and like the cube of the typical slice size, which should 
include spin and Nambu space dimensions in our case. 
Therefore, in a purely single mode wire, all slices have 
size 4, while M mode wires will have slice sizes 4M. 

Once obtained Goo, it can be plugged into the G in 
the Lippmann-Schwinger equation, since the coupling V 
connects the leads only to slice zero, by definition. The 



Multimode, smooth pinch-off 
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FIG. S2. dl/dV in s. N'^S'^S junction for various lengths of 
the superconducting depleted region, Lgd^ and the normal 
depleted region, L^a^ and weak spin-orbit coupling, Lso = 
5La. Different columns feature increasing values of Lgd from 
left to right, whereas different rows feature increasing length 
L^d from top to bottom. 



FIG. SI. Thermal smearing (top) at T = 70mK of the dl/dV 
response in Fig. 4(d), with a section at fixed B — 0.5Bc for 
increasing temperatures (bottom). Note the suppression of 
the ZBA amplitude far below the quantized zero temperature 
value 2e^//i. 



scattering matrix can then be read off by decomposing 
l^a) into right-propagating and left-propagating modes 
in each lead. In our simple case, this can be done by com- 
puting the bandstructure of the lead, but in a more gen- 
eral case the decomposition may be performed by analyz- 
ing g. Additionally, to compute the differential conduc- 
tance dl/dV in an JV channel Normal-Superconductor 
junction 



dI{V) 
dV 



e 
~h 



J\f - TrO 



.t 



- Tr(r^^re/,) 



e=y 



a decomposition into the electron (e) and hole {h) sectors 
of the reflection matrix has to be carried out. 



II. GENERALIZATION TO A MULTISUBBAND 
NANOWIRE 

Applying the above recursion technique to multisub- 
band nanowires is relatively straightforward. It is es- 
sentially achieved by considering a finite width in the 
transverse {y) direction of the wire. The y dimension is 
then discretized much like the x direction, and included 



into the mesh to solve recursively. Some physical de- 
tails must be observed, however. The spin-orbit cou- 
pling in this case should include a dy component. For 
Rashba coupling, for example, this introduces an addi- 
tional j dx ip^{x){—iaaxdy)ilj{x) term in H. Addition- 
ally, one must consider the energy shift of each of the 
resulting bands in the leads due to confinement energy 

J dx tlj^{x) ( ^^ J tlj{x). The lead integration and recur- 
sion procedure is carried out in the same fashion, albeit 
with larger sized slices and more costly matrix inversions. 



III. EFFECT OF FINITE TEMPERATURE 

Finite temperatures can be incorporated into the sim- 
ulation rather easily, since its effect on dI/dV\j^ ^ 
amounts to a convolution of the zero temperature result 
dI/dV\Q and the derivative of the Fermi-Dirac distribu- 
tion: 



dl 
dV 



kBT 



dV 



dl 
dV 



d 



-1 



dVl^e^/^B^' 



(SI) 



In the experiment^, a typical temperature is T ~ 
70mK, which is a rather small scale in units of the super- 
conducting gap A (= 250/iel/), KbT/A ^ 0.024. How- 
ever, as seen in Fig. SI, such smearing has a strong effect 
in reducing the height of the narrow resonances at zero 
temperature. The width of the resonances in units of 
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FIG. S3. dl/dV in a N'^S'^S junction for various lengths Lgd 
and Lj^d and intermedium spin-orbit coupling, Lso = ILa- 



FIG. S5. dl/dV in a NS^S junction for various lengths L^d 
and Ljsid and weak spin-orbit coupling, Lso = 5La. 
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FIG. S4. dl/dV in a N'^S'^S junction for various lengths Lgd 
and Ljyd and very strong spin-orbit coupling, Lso = 0.2La. 




FIG. S6. dl/dV in a NS^S junction for various lengths Lgd 
and LjYd and intermedium spin-orbit coupling, Lso = ILa- 



ksT gives a sense of the thermal suppression of the sig- 
nal. Since, moreover, in multimode wires the width of the 
zero-bias anomaly (ZBA) coming from Majorana modes 
is strongly suppressed by the pinch-off of deeper modes, 
the overall ZBA amplitude at finite temperatures appears 
to be much smaller than the nominal zero-temperature 
2eyh. 



IV. EXPLORING THE PARAMETER SPACE IN 
THE dl/dV RESPONSE 

The different characteristics exhibited in transport 
spectroscopy for spin-orbit lengths of 200nm, as corre- 
sponds to InSb wires, are presented in the main text 
in Fig. 2. Here we extend this analysis (also for spa- 
tially abrupt potentials) by considering wires with weaker 
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FIG. S8. Zero temperature dl/dV at constant Zeeman field B 
below Be, B = 0.5Bc in (a) and B = 0.7Bc in (b). The value 
of Be = y^(/i — t/d)^ + A^ varies as the depletion potential 
Ud is swept, and a transition occurs at 5^ = B, dashed white 
line. When B > B^, the depleted S^ region enters a topolog- 
ical superconducting phase, and Majorana bound states are 
formed. The residual splitting is much smaller than the res- 
onance spacing of the Majorana mode precursors visible at 
B < Bf. 



FIG. S7. dl/dV in a NS'^S junction for various lengths L^d 
and Lp^d and very strong spin-orbit coupling, Lso = 0.2La- 



and stronger spin-orbit coupling. Figs. S2-S4 present 
the results for N^S^S wires, while Figs. S5-S7 corre- 
spond to NS^S wires. Note that, in all these panels, 
there is a certain asymmetry between the -\-V and the 
—V behaviour of the dl/dV response^. This is due to 
the asymmetric electron-hole quasiparticle transport be- 
haviour that is present in systems where the Andreev 
approximation doesn't hold, i.e., in systems in which A 
is not much smaller than /i. This is the case of the kind 
of semiconducting nanowires that we study here, where 
we use /i = 4 A in all these plots (other parameters are: 
Ud = 3.25A, (5 = as corresponds to sharp potential 
profiles, and Up ~ 30-70 A). Note that this effect is also 
present in the dl/dV calculations of the main text. 

The following general features are apparent as spin 
orbit strength is increased (spin-orbit length Lso is re- 
duced): (i) the spacing and the width of Andreev bound 
states (ABSs) increase; (ii) the band edge at large mo- 
mentum (green lines) becomes less sensitive to Zeeman 
coupling B (spin becomes more polarized in plane by 
the spin-orbit coupling); (iii) Although the formation of 
Majorana modes and their corresponding ZBAs are not 
affected by the strength of the spin-orbit coupling for 
B > Bc^ the MBS pairs forming in the depleted region 
for Bf < B < Be exhibit a splitting due to wavefunction 
overlap that grows with spin-orbit coupling strength. On 
the other hand, a weak spin-orbit coupling has the dis- 
advantage that the excitation gap above the Majorana 
modes is diminished, which has important implications 
for their coherent manipulation for quantum computa- 
tion purposes. If Lso is smaller than La, and the de- 
pleted region is long enough so as to suppress the over- 
lap splitting, one could achieve quasidegenerate Majo- 



rana pairs isolated from the quasiparticle continuum by 
a gap equal to a full A, improving the chances to per- 
form coherent non-Abelian braiding transformations of 
the wire's ground state at Bf < B < Be. 

Note once more that zero energy ABS anticrossings 
may coincide with Majorana pairs in the depleted region 
(above Bf) in the NS^S junction case, bottom rows of 
Figs. S5-S7. Such circumstance is not possible in the 
N'^S'^S junctions of Figs. S2-S4. For long enough Lgd 
lengths (bottom-right panels), the MBSs around zero en- 
ergy decouple from the pinch-off probe, and the ZBA 
signal is suppressed. This wavefunction repulsion effect 
between ABSs and MBSs as seen in the experiments, 
only arises between resonances belonging to the same 
subband, at least in the absence of subband mixing ef- 
fects such as disorder. 



dl/dV RESPONSE UNDER A SWEEP OF 
THE DEPLETION GATE 



Experiment in Ref. 2 includes some results on the 
dl/dV response in NS^S junctions as the depletion gate 
Vd is swept at constant Zeeman field B (see e.g. Fig. S5 
in their supplementary material, Vd is called Vi there). 
The traces revealed what appears to be a split Majorana 
ZBA at small values of Vd- The authors point out this 
splitting could be related to overlap of Majorana bound 
states, but that more work needs to be done to confirm 
this interpretation. Our simulation, see Fig. S8, shows 
that the type of splittings obtained in the experiment 
are probably not due to Majorana overlap, but rather 
to a transition from B^ > B (to the left of the dashed 
white line) to B^ < B (to the right of the dashed white 
line) as Ud is swept. Hence, the split ZBA is actually 
the remnants of the gap closing that is to be expected 



le 



in the standard picture of Majorana state formation, in 
the form of a low-lying ABS that is the precursor of a 
Majorana pair in the depleted zone. The actual splitting 
due to Majorana overlap occurs when the split resonance 



merges into a much narrower ZBA, whose residual split- 
ting cannot be resolved experimentally at the working 
temperatures. 
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